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ABSTRACT 


This paper presents the numerical calculation of the Wannier functions (a,(z)) of 
the GaAs/Alo.2sGao7sAs superlattice structure. The Wannier functions are linear 


combinations of the Bloch functions (y,, (z)) that can be viewed as a convenient 


mathematical instrument to get around the lack of orthogonality of the tight binding 
formulation and are useful tool when the position of an electron has physical importance. 
However, except for finding a variational principal for the energy levels of one 
dimensional crystals in terms of The Wannier functions and simple cubic lattice 
calculations for instructional reasons, models have not been calculated or plotted showing 
the theory can be expanded to more complex superlattice structures. We develop an | 
algorithm to calculate numerically and plot the Wannier functions of the 
GaAs/Alp 2sGao7sAs superlattice structure and hence prove that these functions can be 
calculated even for complex structures. By using the plots of our numerical modelling, 
we displayed the peculiar properties of the Wannier functions: 

a. Wannier functions are real. 

b. They fall off exponentially. 


c. They are either symmetric or anti-symmetric about z=0. 
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I. INTRODUCTION 


Semiconductor superlattices are a novel class of man-made solids which consist 
of a repetition of alternating layers of two or more different semiconductor materials. By 
adding artificial periodicity to the multi-layered structure [Ref. 1], the electronic nature of 
the superlattice can be controlled through the constituent material’ s crystalline 
orientation, layer thickness and material composition. The result of this periodic layering 
process creates semiconductor superlattices which present exciting optical and electronic 
characteristics which cannot be found in existing semiconductor crystals [Ref. 2]. 

From the beginning of the study of this new class of solids, most of the attention 
and research has been concentrated on new generations of optoelectronic device 
applications. One of the promising device implementations is fast and highly sensitive 
InAlGaAs/InAlAs superlattice avalanche photodiodes for use in 1.3 to 1.6 um optical 
communications [Ref. 3]. Other areas of research include the study of superlattices for 
the future generation of far infrared detector imaging arrays with wavelength A>12 um 
[Ref. 4] and the study of high optical confinement effects [Ref. 5] leading to 
investigations of using superlattices to reduce the threshold current density and enhance 
the device stability of semiconductor lasers. 

Even though successful device implementations have been reported, and a world- 
wide research effort is underway in this new field of study, many problems remain 
unsolved. These problems include everything from the preparation of the superlattices to 
understanding the transport and optical phenomena introduced by these novel materials. 

In this thesis, the Wannier functions (a, (z-d;)) of a GaAs/Gap7sAlo.2sAs 


superlattice is numerically calculated and plotted with the new lattice constant in the 





growth direction z set by dgaastdaicass (Fig. 1). The Wannier functions can be viewed as 


V(z), eV 
Superlattice 
unit cell 
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Figure 1: Conduction band edge profile of GaAs/Alp.7sGao2sAs superlattice with the new 
lattice constant d = daicaas + dcaas Where daicaas and dgass are are thickness of the 
corresponding semiconductors. 


a sauvenicat mathematical instrument to get around the lack of orthogonality of the tight- 
binding formulation . It is also a useful tool in understanding the spatially localized 
nature of electrons when they are near an impurity state. Except for theoretical studies on 
finding a variational principal for the energy levels of one dimensional crystals in terms 
of the Wannier functions [Ref. 6], and simple cubic lattice calculations for instructional 
reasons [Ref. 7], models have not yet been calculated or plotted showing the theory can 
be expanded to more complex superlattice designs. This thesis will numerically calculate 
and plot the Wannier functions of GaAs/Alo.25 Gao.7s AS superlattice structure to prove the 


theory beyond the simple cubic problem and will display the properties of these 


functions graphically. 








i. THE WANN IER FUNCTIONS AND THEIR PROPERTIES 


In a one dimensional crystal, the mathematical basis for the motion of an electron 


can be expressed by the Schrodinger equation as 
hh? 2 
- 2S aveto yay (1) 
m 


where V (z) is a periodic potential function of the crystal in the growth direction z with 
the period of lattice constant c, ie., V/z)=V(z+c) and m is the mass of the electron. Bloch 
has shown that the eigenfunctions of the one-electron Schrodinger equation, or Bloch 
states, for this periodic potential (Eqn. 1) have the mathematical form [Ref. 8] 
Wy (Z) =e Uy (Z) (2) 

where u,,(z) is a function with the periodicity of the potential. The subscript n stands for 
the n™ energy band to which the eigenfunctions belong and the subscript k is the wave 
vector in the growth direction z. 

In the method of the tight binding formulation, using atomic wave fanctions as a 
basis set has a practical advantage over Bloch states in qualitative calculations [Ref. 9]. 
The actual calculations with those functions, however, create difficulties due to the 
atomic wave functions being centered on different atomic sites that are not orthogonal to 
each other. Furthermore, they are not necessarily orthogonal for different k values either. 

Wannier has introduced a set of localized wave functions, a,(z-c;), which are 


linear combinations of all the Bloch states y,, (z) for a given energy band, and have the 


advantages of orthogonality [Ref. 10]: 


a,(2-6)) == De Wy (z) (3) 





where 
1 
a,, (Z) =p 2m (Z) (4) 


N is the number of lattice sites in the one dimensional crystal and subscript j stands for 
the lattice site of which we are going to calculate the Wannier function. The summation 
in Equation 4 extends over all values of the first Brillouin zone (BZ). The Bloch function 
y.,(z) is continuous in k, and is the same function of & at k+K where K is a reciprocal 
lattice vector, i.e., the Bloch functions are periodic in the reciprocal lattice [Ref. 11]. 
Since the Wannier functions are functions of distance from the direct lattice site j, we 
may think of Equation 3 and 4, i.e., yw, (z) and a,(z-c)) as the Fourier transform of each 
other. The Fourier-transform of Equation 3 defining the Wannier functions, can now be 


calculated to obtain the Bloch functions: 
] ike ; 
(Z)=—= ) e /a,(zZ-C¢;) (5) 
Wax iy 2 | 


The Wannier functions can show orthonormality properties in contrast to the tight 
binding formulation, because the Wannier functions centered on different lattice site 


and/or or belonging to different energy bands are orthogonal to each other : 


fa, (z-¢))a, (z-¢;)dz = wee [Va OV, 4: dz 


1 i(ke;—-k ¢;) 
ga? as Sim 


= Om cc (6) 











where we have invoked the summation relation 


y etre — NO... (7) 


k 
and the orthogonality of the Bloch functions. Hence, the Wannier functions form a 


complete orthonormal set, and the completeness can be shown by summing over all 


energy bands and all lattice sites, 


SY) a, (z-¢;)a,(z —c;) =6(z-z). (8) 
n,j 


For a one dimensional periodic potential with no degeneracies of energy bands, It 
has been rigorously shown that the Wannier functions will decay exponentially at large 
distances [Ref. 12]. For our purposes, we can easily show this localized nature of the 
Wannier functions qualitatively. If we replace z with c; in Equation 3 where c¢; is a lattice 
vector, then we get 


—ik (c;-¢;) 


a,(¢\-¢) = De v,,(2) (9) 
k 


When | c; — c; | becomes large the exponential part of the summation oscillates rapidly and 
contributes vanishingly to the summation. Hence, we can expect the Wannier functions 
of one energy band to fall off exponentially at large distances and to be identical 
functions which are centered at ite of the lattice points. 

One of the most interesting properties of the Wannier functions is their reality. 
Regardless of the phases of the Bloch functions, Wannier functions are real. W. Kohn | 
has discussed the reality problem of the Wannier functions in detail [Ref.12]. We are not 
going to prove this peculiar aspect of the Wannier functions in this paper, but we refer the 


reader to W. Kohn’s paper. 


With the aforementioned properties, the Wannier encore can be considered a 
useful mathematical tool when the position of an electron has physical importance, e.g., 
in the Heciaiton of the motion of an electron around an impurity site. Other important 
areas of application include deriving a transport theory for Bloch electrons [Ref. 13]. 
Also, using Wannier functions in the calculation of energy bands has been discussed but 


no practical application emerged [Ref. 14]. 








Ill. SUPERLATTICE EFFECTIVE-MASS SCHRODINGER 
EQUATION 


Before presentin g our method of numerical calculation of the Wannier functions, 
we are going to derive the Schrodinger equation for the GaAs/Alo75Gao.25As superlatice 
structure (Fig. 1) by using the effective mass theory with the Kronig-Penney model. This 


approach describes the electron with an effective mass m, where subscript i corresponds 


to the constituent layers of the superlattice, and approximates the superimposed 
perturbing potential (Superlattice bandedge profile) as periodic square wells (Fig. 1). 


If H is the Hamiltonian for an electron in a one dimensional perfect crystal, then 


Ay, (2) = E,(k)Wy. (2) (10) 
where 
hi 0? 
H = “on : (11) 


and E(k) are the eigenvalues of the Bloch functions y,,(z), namely the band structure 
function. Since band structure function E(k) is a non-degenerate and continuous 


function of k with the periodicity of the crystal in the reciprocal space (k-space) [Ref. 15], 
we have the possibility of defining the Fourier expansion of it as we did for the Bloch 
functions (Eqn. 5): 


E,(k)= Ee 2 (12) 


where the summation is over the lattice sites of the host material in the growth direction 


(z), and EF, (c)is the Fourier expansion coefficient. 


Since the band structure function E,(k) is an analytic function of k and the wave 


functions y,,(z) have the periodicity of the crystal, we can assume the replacement of 


the perfect crystal Hamiltonion (H) by the operator E,,(—7 a which follows from 
Z 
E_(k) when Hermitian differential operator — Pe is used in place of k [Ref. 16] where i 
Z 


is ¥—1. We will show that H and E, (-i = have the same eigenvalues and 
Z 


eigenfunctions. From Equation 10, we have 


2 
E, Ci. (2) = VE, (e ®Vy (2) 
0 1,0. 
= Dbl res +565. + Wal 


1,0 
= DE, of Ya iZ)-7 CW, (Z) +50 aU (Z)+ -| | 


The term on the right hand side in the brackets is nothing but the Taylor expansion of the 


Bloch function y,,(z +c). Therefore, we obtain 
i) 
EAS Wan (2) = LEW ulZ +0) 
= YE, (c)ey_, (Z) 


=E,( Wy (Z) - (13) 
We have used the Fourier transform of the band structure function (Eqn. 12) and Bloch’s 
theorem [Ref. 17]: 


Wy (Zte)=ehWy (Zz) (14) 











As can be noted, the periodic potential of the crystal structure V(z) in Equation 11 is not 
explicitly present in this formulation (Eqn. 13) where it is incorporated with the usual 
kinetic energy operator, 


h? 0? 


hog 15 
2m dz” “) 


The unperturbed Hamiltonian defined by the modified kinetic energy operator 


2 2 
E(-i S would be far more complicated than the kinetic energy operator ne Q 
iz , 


= 52 , if 
we tried to include the full k-dependence of the band structure function E,,(k) [Ref. 18]. 
However, many aspects of the electronic properties of semiconductors can be obtained 
from a small number of excited electrons found at the bottom of the conduction band and 
holes created in the absence of these electrons at the top of the valence band. Hence, what 
we need is the form of the band structure function for these carriers over a small region of 
k-space around the band extrema. This implies that the bandstructure function can 
generally be drawn as a quadratic function of k. Besides, for most of the semiconductors, 
E_(k) is approximately parabolic near an extremum, as for the conduction band edge of 
GaAs and Alo7sGao2sAs [Ref. 19]. For the present problem, by restricting our treatment 
to the band extremum E, occuring around the wave vector value k =k, =0 which is the 
case for direct gap semiconductors like GaAs, we take the Taylor series expansion of 


E(k): 


E(k) = Ey + Bk =k) +5 By (K — Ky)? # oveees (16) 


where the expansion coefficients are 





aE 
Baal ¢ (17) 
Ok’ poe 

the subscript i representing the order of expansion. Since the Taylor expansion of E-k 


_ relation around the band extremum is relevant to the problem, the first order term 


B,(k—k,) will vanish by construction, and higher order tems can be neglected. For 
dimensional reasons, £, must have the units of h’ /m. Therefore, by defining an 
effective mass (m") of the electron and replacing the curvature of the E-k relation with 


hi? /2m*, the band structure function near the extremum point can be characterized as 


2 





E,(k) = Fy+— kes. , (18) 


m 
Once again, by replacing the canonical substitution — Lae in place of k in Equation 18, 
| Z 


we obtain the effective mass Schrodinger equation for a bulk crystal, 


Fare i Bs al) = E,(K)W x (Z) > (19) 
where the effect of the eneeal potential V(z) is absorbed in the parameter m 

Now that the Schrodinger equation for the perfect crystal has been formulated, we 
can introduce our perturbing potential U,(z) shaped by the GaAs/Alp 25Gao75As 
superlattice in the growth direction z. By applying the Kronig-Penney model depicted 
earlier, the superimposed superlattice potential can be expressed as 


U. (z) _ | ia - Em = AE. lz —nd))d gaa, /2 


=0 le—nd\deu,/2 (20) 








where ES“ and E4°* are the bottom of the senltitiice bands of GaAs and 
Alp 2sGao75As respectively and d is the superlattice constant defined earlier. Equation 20 
explicitly implies that U,(z) is a periodic function of z with the periodicity d, i.e., 
U (z)=U,(zt+nd) . (21) 
By defining the perturbing superlattice periodic potential U,(z) , which is 
superposed over the intrinsic crystal potential of the host material, the Schrodinger 


equation for GaAs/Alo.25Gao7sAs superlattice can be then written as 


Z 2 
-A Sw, We (Z)=E, (KW (2), (22) 


where we moved E, to the right-hand side of the Equation 22 and defined €,(k) which 
is the difference between bandstructure function E_(k) and band extremum E,. This 
rearranging of energy terms will give the energy values from the edge of the conduction 
band of the constituent layer (Since EC“ < E“°S , in this context it is GaAs) which is 
_approximated as a quantum well by the Kronig-Penney model. The spatial dependence on 
the growth dnection z is brought in to the effective mass (m'(z)) since the materials 


forming the layers of the superlattice have different values of effective mass. 
Since the Bloch functions are orthogonal for different energy values which are 


real, it requires the hermiticity of the Hamiltonian (Energy operator), 


h> 9d? 
H=- ——+U : 23 
2m’ (z) dz? sa on 





But, this form of the Hamiltonian is not accurate when the effective mass has spatial 


e « : e re] a 2 
dependence , since the canonical operator — i does not commute with z, which 
Z 


11 


terms. For example, let us consider an 





0 
requires the correct ordering of —i— and 


dz m (Z) 
interface between the two constituent layers (GaAs and Alo.25Gao.75As), Say at Zp). The 
wave function y,,(z) has to obey the usual boundary conditions (Ref. 20] to guarantee 


current conservation across the junction, 


Vu (Zp = WarlZo (24) 


OW (Z) _ Q Wy, (Z) (25) 
dz Gas dz ZAlGahs 








where 2°“ and 24'S" are the sides of the junction in GaAs and Alo.2sGao7sAs layers 
respectively. On the other hand, the layers have different effective masses [Ref. 21}. So, 
the boundary condition for matching the derivatives of the wave functions (Egn. 25) must 


take the effective masses into account: 


1 ow,(z)) — 1 OW,(z) (26) 





* 0 i, * ad 
Mads Z yA m AlGaAs x Z ans 


where m,,,, and M Gass ate effective masses of the electron in the corresponding 
semiconductors. Since the derivative in Equation 26 is apparently the Hermitian 
differential operator, the hermiticity here also guarantees the current conservation. Asa 
result of these consideration, the correct form of the Schrodinger equation with a spatially 


dependent effective mass is given by 


] 
72 : OV) 4 7 (zyy, (2) =E (KW y(Z), = (27) 
Ozm(z) az 





which assures the hermiticity and guarantees the conservation of current at the interface. 


12 











IV. THE WANNIER FUNCTIONS IN FINITE DIFFERENCE FORM 

In this chapter we will present the numerical approach which is utilized to obtain 
the Wannier functions of the GaAs/Alp2sGag7sAs superlattice structure. The algorithm is 
then built in two steps: The first step is to solve Equation 27 for the wave functions 
y_,(z)and the eigenvalues €,(k) (Bandstructure function). In the second step we ni 
need to use Equation 4 to calculate and plot the Wannier functions; i.e., for each band 
index n, we sum the wave functions over all values of k. 

We are going to use finite difference approximation [Ref. 22] to solve the | 
superlattice effective mass equation (Eqn. 27) for the wave functions. In this method, we 
evaluate the values of y/,,,(z) at discrete points z; which are obtained by dividing the 
superlattice unit cell into finite number, N-1, of intervals 4z defined by 

AZ=Zj-2, (28) 
The first derivative in the finite difference form at the point z; may be written as 


oO 
Oz 


(z, + Az)—W,, (z; — Az) 
Vn (Z)|,, ee (29) 


The expression for the second derivative can be obtained by applying Equation 29 


successively [Ref. 23]: 












a (z)}) = 
dz m'(z) dz" ™ |, 
] 1 1 ] 
ee a ee (2, ) 4 Ve) 
- Ce j+2 Mm (Zig) (Z4) eae dee Ca j-2 





30 
(2Az)? (30) 


Before writing the eigenvalue matrix equation for Equation 27, we have to take 


care of the boundary conditions along the growth direction z. In the artificial periodicity d 


13 


of the superlattice, the wave function must satisfy the Bloch theorem (Eqn. 14) [Ref. 17]: 
yw (ztd)=e"y,, (z) : (14) 
In the finite difference matrix equation, we replace w,,(Z)) and Y,,(Zy4,) according to 


Equation 14. Since we have expressions for each mesh point z; in the superlattice unit 
cell, Equation 27 can be rewritten as an eigenvalue equation. By setting 


2Az = @=dI(N —1) and an replacing the constants in Equation 27 with 





T hn Tr T hn” i 
AlGaAs — 2m’ eure” ~~ "A GaAs ~ 2m 4.0" — “G 
nk? nk? 
Vaicaas = AE, +5 : =a Vous =5>= — = Ve (31) _ 
MM aiGaAs Maas 


and using Equation 30, we get the eigenvalue matrix equation for each value of k: 


—2T,+V, T, 0 O 0 0 0 O T,e“ 
r, —-27,+V, IT, 0 0 0. 0 0 0 

0 lr, 0 0. 0 0 0 

0 0 ae iE 0 0 0 0 

0 0 0 IT, -21,+V, I 0 0 0 

0 0 0 0 hi; -~2T,+V, Tz 0 0 

0 0 0 0 0 A; 0 

0 0 0 0 0 0. r, 

je 0 0 O 0 0 0 YT, -27,+V, 
Wx (Z) Wi (Z) 
Wn (22) | W,1(Z) 
=é,(k)) . (32) 
Wx (Zy) Vin (Zy) 
14 








Once we have obtained the wave functions for each value of the reciprocal wave 
vector k for a specific band index n, we will create a subroutine to sum the wave 
functions over all the k values in the first Brillouine zone (BZ) (Eqn. 4) to calculate and 


plot the Wannier functions of our superlattice structure. 


15 
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V. RESULTS AND CONCLUSIONS 

The computer codes in the Appendix A afd B have been developed for the 
solution of the effective mass Schrodinger equation for the GaAs/Alp2sGag75As 
superlattice consisting of, alternately, Alo 25Gao7sAs and GaAs layers. For the purpose of 
this paper, we have chosen the thickness of the Alo.2sGao7sAs layer as 6 nm and because 
of our familiarity with it, the Matlab programming language has been used. 
A. BANDSTRUCTURE FUNCTION 

Before presenting the objective of this paper, i.e., plotting the Wannier functions 
of the GaAs/Alp.25sGao.7sAs superlattice, we wanted to make sure that the eigenvalue 
matrix equation (Eqn. 32) that we formed as the discretization of the superlattice 
Schrodinger equation (Eqn. 27) is accurate. The easiest and the most effective way to do 
so would be to compare the eigenvalues of the matrix equation (Band structure function) 
to another method (which has already proven itself to be correct) of finding the band 
structure function of a superlattice structure. 

We first created a subroutine in Appendix A for the evaluation of the eigenvalues 


é,(k) of the matrix equation (Eqn. 32) which gave us the bandstructure function of the 


superlattice. Next in the second subroutine, we used G. Bastard’s analytic superlattice 


dispersion relation equation [Ref. 24] 


Cos(kd) = Cos(k,d, )Cos(k,d,)- : + FE bite )Sin(k ,d, ) (33) 


ke he 


17 


2€, (k) Means 28, (KM car 
ke = Keaas oor ae » hy =Kaicaas = JAE. ae maa 
dg — A caAs : d, = A stGaAs (34) 


to compare with the eigenvalues of the matrix equation. Equation 33 and 34 are implicit 
equations relating the energy values €, (k) to the wave vector k. 

Figure 2 shows the bandstructure function obtained by evaluating the eigenvalues 
of the matrix equation (Eqn. 32), while Figure 3 is the solution of the Equation 34 
satisfying the dispersion relation (Eqn. 33). By inspection, it is clearly seen that the two 
bandstructure figures are in good agreement with each other. That is to say, the 
eigenvalue matrix equation formed as the discretization of the superlattice effective-mass 


Schrodinger equation (Eqn. 27) is accurate. 


BANDSTRUCTURE OF GaAs/AlGaAs(Al%25) SUPERLATTICE 
0.45 


AE. = 0.2026375eV 


d cass = OnM 


oe 4 aiGaas = 3M 


09 YO ST by 
2) 
N 


eV 





3 -2 1 O 1 2 3 

k (First BZ), I/m 
Figure 2: Bandstructure function ¢€,(k) obtained by evaluating the 
eigenvalues of the matrix equation (Eqn. 32) where AE, = E“°“ — E°™s, 
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BANDSTRUCTURE OF GaAs/AlGaAs SUPERLATTICE (Eqn. 33) 


0.45 
0.4 AE, = 0.2026375eV 

" 0 35 d Gass = 6nm 
e , d siGass = 3M 
r 
g 0.3 
y 

0.25 
eV 





x 10 


k (First BZ), I/m 


Figure 3: Bandstructure function €,(k) obtained by solving the G. Bastard’s _ 
superlattice dispersion relation equation (Eqn. 33 and 34) where 
AE. = EAGaAs _ Gods” 


B. WANNIER FUNCTIONS 


After confirming the accuracy of the numerical modeling of the superlattice 
Schrodinger equation (by comparing with the G. Bastard’s superlattice dispersion 
relation), the computer code in Appendix B implements an algorithm which calculates 


and plots the Wannier functions of the GaAs/Alo 2sGao7sAs superlattice structure. At 
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this phase of our numerical modelling, we deliberately skip the normalization of the 
Bloch waves. There are three reasons for this: First of all, the Bloch functions can easily 
be calculated and have been obtained by solving the Schrodinger equation for each value 
of k and n [Ref. 25]. Secondly, we a not concerned with the normalization of the Bloch 
functions, since having arbitrary amplitudes for the Wannier functions of any band index 
(n) is immaterial for the purpose of this paper, i.e., plotting and pointing out the 
properties of the Wannier functions. The third reason which is the limitation of our 
modelling is that normalizing each Bloch function of every superlattice unit cell would be 
formidable to carry out numerically if you discretized 300 mesh points for each 
superlattice site (z;) and 100 mesh points in the first BZ (x) which brings out 30,000 
Bloch functions just for the first band index (n=1) to deal with. We are going to point out 
this computational difficulty in the conclusion ia final remarks subchapter. 

In Chapter II, we have with inclusive detail proved that for each band index (7) 
there exists only one Wannier function. Since for most of the optical and electronic 
applications the first two bands of the bandstructure function are utilized, in the computer 
code (Appendix B) we numerically calculated the Wannier functions of those bands, i.e., 
n=1 and n=2. Favorably, plotting the Wannier functions for these two bands suffices for 
the purpose of this paper. Furthermore, we developed the computer code in such a way 
that by changing the band index entry in the program we can obtain the Wannier 
functions of any band index that we are interested in. 

Figure 4 shows the Wannier function of the first band (a;(z)) extending to the 
right neighboring superlattice site. By inspection, it is clearly recognized that the 


Wannier function (a;(z)) is localized and centered at the superlattice site which is shown 
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WANNIER FUNCTION OF GaAs/AiGaAs(Al%25) SUPERLATTICE 


Conduction band edge profile of superlattice 


daiGaAs= 3 nm 


dcaAs =6nm 


Amplitude (unitless) 





O 0.5 1 1.5 2 2.5 3 
Z, (nm) 


Figure 4 : The Wannier function of GaAs/Al»-; Gan75 As superlattice for the first 
band index (@;(z)). 


with a thick line in Figure 4. The asymptotic behavior is also observed easily as the 


Wannier function decays down in the neighboring suverlattice unit cell since as lz —d ; 


increases (Eqn. 9), the Bloch functions destroy each other by interference which we had 
proved thoroughly in Chapter II | See the discussion after Equation 9). 

We shall now show that the Wannier functions are real. That is to say, if we tried 
to plot the ‘faeenaarput of the numerically calculated Wannier function of the first band 


(or any band), we should get an empty figure. However, since the numerical approach 
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(finite difference method) that we used to solve the superlattice Schrodinger equation 


(Eqn. 27) is just an approximation, we would expect some amplitude as error. Fi gure 5 


x 10> WANNIER FUNCTION (Imaginary) OF GaAs/AlGaAs(Al%25) SL 


Amplitude (unitless) 





1.5 2 2.5 3 


0 0.5 1 
_ Z, (nm) 


Figure 5: Imaginary part of the numerically calculated Wannier function of 
GaAs/Alo25 Gao.75 As superlattice for the first band index (n=1). 


shows the imaginary part of the Wannier function of the first band (a;(z)). Indeed, if we 
note the order of amplitude which is around 10°’, we can confidently state that the 


Wannier functions for theGaAs/Alp 25Gag.75As superlattice is real. 
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Figure 6 depicts the Wannier function of the superlattice (SL) structure for the 


band index n=2 (a2(z)). The localised nature and the asymptotic behavior are obviously 


recognized as we did for the a;(z). Again, the Wannier function a2(z) is centered at the 


WANNIER FUNCTION OF GaAs/AlGaAs(Al%25) SUPERLATTICE 


Amplitude (unitless) 


Conduction band edge profile of superlattice 


daicaAs= 3 nm 


doaas =6nm 





O 0.5 1 : 1.5 2 2.5 3 
zZ, (nm) 


Figure 6: The Wannier function of GaAs/Alo.2s Gag75 As superlattice for the 
second band index (a2(z)). 


same superlattice unit cell as we have verified for a;(z) in Figure 4. The important aspect 
of Figure 6 is the anti-symmetry of the Wannier function for the band index 2 (a2(z)) 
about z=0 where we can easily see the symmetry of the Wannier function for the band 
index n=1(a;(z)) in Figure 4. W. Kohn has given a comprehensive discussion about the 
symmetry properties of the Bloch and the Wannier functions [Ref. 25]. Basically the 


(anti-) symmetry property of the Wannier functions about z=0 can be accounted for by 
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“the parity of the Bloch functions which could be easily perceived if we revisited the basic 


quantum mechanics problem of square well. 


The reality of the Wannier functions can again be shown if we plot the imaginary 


part of a2(z) which is resulted from our numerical approach (Fig. 7). Note the 


x 10° WANNIER FUNCTION (Imaginary) OF GaAs/AlGaAs(Al%25) SL 


Amplitude (unitless) 





2.5 


are) 0.5 1 1.5 2 
Z, (nm) 


Figure 7: Imaginary part of the numerically calculated Wannier function of 
GaAs/Alo25 Gao75 As superlattice for the first band index (n=2). 


amplitude of the numerically calculated imaginary part of a2(z) which is around 10°). 


C. CONCLUSION AND FINAL REMARKS 
In this thesis, we have developed an algorithm to numerically calculate the 


Wannier functions of a GaAs/Alo 25 Gao.75 As superlattice structure. The algorithm has 
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been built in two seperate steps : First we wanted to show the credibility of our 
numerical modelling of the superlattice effective-mass Schrodinger equation (Eqn. 27). 
We did that by comparing the bandstructure function modelled by our method with the 
one obtained by using the G. Bastard’s superlattice dispersion relation (Appendix A). 
We proved that our modelling with the finite difference form is a very good 
approximation for the discretization of the superlattice effective-mass Schrodinger 
equation (Figures 2 and 3). In the next step, after sonanmine that our numerical 
modelling is reliable, we created a routine (Appendix B) for solving the eigenvalue 
matrix equation (Eqn. 32) for the Wannier functions of the GaAs/Alo.25 Gao7sAs. 

The results we obtained are in good agreement with the theory of the Wannier 
functions which we have discussed in detail in Chapter I. Although, in the past there 
were no attempts to calculate these functions for bulk crystals, we have demonstrated 
that the Wannier functions can be calculated even for a complex superlattice structure. 
By plotting the numerical simulation of the Wannier functions for the first couple of 
allowed bands ( n=1 and n=2) of the superlattice structure , we displayed the peculiar 
properties of these functions: 

a. Wannier functions are real (Figures 5 and 7). 

b. They fall off exponentially (Figures 4 and 6). 

c. They are either symmetric or anti-symmetric about z=0 (Figures 4 and 6). 

Even though the results we reach conform with the theory, the numerical 
approach we have utilized create two shortcomings that we have to deal with. The first 
one which we pointed out before is the normalization problem of the Bloch functions. 


Since for every band index n and for each value of k there is a Bloch state (y,, (z)), 
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normalization for a continuous k values is a formidable task to do, e. g., just for the first 
band, the computer code in Appendix B has simulated 16800 Bloch functions for discrete 
values of k. Besides, we have to note that we just extend the Wannier functions to the 
next two superlattice sites to show the asymptotic behavior. But, as we have stated 
before, a purpose of this paper it is unnecessary to normalize the Bloch states. The 
second problem we encounter is about the effective masses. In order to show all the 
aspects of the Wannier functions we have had to use a uniform effective mass throughout 
the simulation distance. The reason for that is the same as the first one: We could not 
discretize the superlattice site small enough to show the kink at each interface the 
Wannier functions should have had that would be due to the different effective masses at 
the interface. We used the smallest mesh distance that we could, but all we obtained was 
discontinuities at the interfaces. 

Finally, we want to call attention to some of the future studies that we think will 

support the theory of the Wannier functions: 

a. We have used Bloch states in the eigenvalue matrix equation to obtain the 
bandstructure function of the GaAs/Alo25 Gao.75As superlattice structure. An 
algorithm can be developed to use the Wannier functions in the calculation of 
the energy bandstructure functions of the superlattices. 

b. Al,Ga;.As bulk alloy is a direct-bandgap semiconductor if x<0.45 [Ref. 26] 
There is a crossover from direct to indirect-bandgap structure at about x=0.45. 


The Wannier functions of such a structure can be studied in the future. 
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c. The GaAs/Alo.25 Gao7sAs superlattice structure is a type 1 superlattice. The 
Wannier functions of type 2 and type 3 superlattices can be studied in the 
future. 

d. If Wannier functions of all types of superlattice structures are obtained, alts 
an electronic transport theory can be jeucionel based on the Wannier 


functions. 
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APPENDIX A: COMPUTER CODE IN MATLAB PROGRAMMING 
LANGUAGE FOR NUMERICAL CALCULATION OF THE 
BANDSTRUCTURE FUNCTION 


format long g 


ELEEEEESSLLSEEESEEEEELELESEEEEEEEEELESEEEEEBSEEELEEESELESEEESSEESESESSS 
SELSSEEESETEEEESSEEEESF INPUT PARAMETERS  %%3S8SE3S3SSSSSSEBBSBEESE 
EXSELEESELSELESEESESSEESESSESEESESSESELSESTESESSSSVESESSESTESESSESSES 


EgG=input (‘Energy gap for GaAs (eV)’); 
hh=input (‘barier thickness (nm)'); 
LL=input(’well thickness (mm) ’); 
x=input(’Al concentration (x)‘); 


EIESEREEYEEEEELESSESEEYSSEESETEESEESEELESEEEESGSSEESEBEESS 
SSSEREEEES¢EEE% DEFINING THE CONSTANTS %8S3%3S3%33b3%4S%% 
EESSESSEBLESESEBELESEBEEBEBSEEEESSEESESSETEESESSSSEBSSBSESBE 


hbar=6 .5821le-16; 

m1=0.067; 

me=(0.511e6) /9e16; 

m2= (mi+0.0838*x); 

EQA=EgGt+1.247*x; 

Ec=(0.65* (EgA-EgG)); SCONDUCTION BAND OFFSET 


h=hh/1000000000; 
L=LL/1000000000; 

d=h+L; 

D= (h/2) +L; 

s=3.14/d; 

in=100; 

adz=d/in; 

dg=floor(s/50); 

w=- (hbar*2) / (2*me* (dz*2)); 


ESSESEELEEESEELESESESEEESEEEEEEESEESESESSESSESEESESSESESSES 
EEESSESEELELESEESTESSE HAMILTONIAN (H ) SSESESSSBLSSSSSESS 
EELSLEESEELTEEEESEEELELEELELEEELEEEELEELELEELEVESEEEEEESEES 


for l=l:in 
Zt LjyeHdz*1i-; 
if z(l)>h/2 & z(1)<D 

V(1)=0; 
HH (1)=V(1)-(2/m1) *w; 
m(1)=mi; 
Ww (1) =w/ml1; 

else 
V(1)=Ec; 
HH(1)=V (1) -(2/m2) *w; 
Ww (1) =w/m2; 

end 

end 


Ww (floor (in/6) )=(2*w) / (m1+m2) ;Ww(floor (in*5/6))=(2*w) /(ml+m2) ; 
HH1=diag (HH) ; 


W=Ww (1: (length (Ww) -1)); 
Wl=diag(W,1); 
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W2=diag(W,-1); 
HH2= (HH1+W1+W2); 


252593999 22 F SILI E LIES EEELESEEELESEEBESEESBESESEESES 
GILLES ESSS FIRST BROULINE ZONE (BZ) %%3%%3%%3%3% 
LALLA LIIIIRL ELLY LEY LS SSERSESEESLELESEESLESEEEESEESSSES 


for k=1:50; 


ql (k)= = (k) al © (0 

Q1l=(w/m2) * (cos (qi (k) *d)- _i*sin (ql (k) *d)); 

Q2=(w/m2) * (cos (qi (k) *d)+i*sin(ql(k)*d)); 

Qll=diag(OQ1, (in-1)) ;Q22=diag(Q2,-(in-1)); 


H= (HH2+011+022) ; % HAMILTONIAN WITH BOUNDARY CONDITIONS 


kk=e1ig (H); & BNERGY EIGENVALUES 
ee=sort (kk) ; 


EE1 (k)=ee (1); 
EBE2 (k) =ee (2); 
BE3 (k)=ee(3); 


BAND INDEX n= 
BAND INDEX n 
BAND INDEX n= 


dP of ge 


end | 
qql=-1*flipir (ql); 
q=[fqql ql]; 


EE11=fliplr(EE1) ;E1=({EE11 EE1); 
EBE22=fliplr(EE2); E2=[EE22 EE2]; 
EE33=fliplr (EE3) ;E3=[EE33 EE3]; 


LIRLELIYSLLYISAILIS SI RSI SEALS A ELSERSEE SESS ESSEBESS SESS 
&£2A29E% PLOTTING BANDSTRUCTURE FUNCTION %%% TESEESESESS 
3%%$%%%% WITH EIGENVALUE MATRIX EQUATION SECELESSEREEESE 
LYIBLGLAI~LYALELIIL AAG EEE LERLESELEEERLEEESESEETSEBES 
figure(1),plot(q,E1,’k’,q,E2,'k’,q,E3,'k’); 

axis({-s s 0 .45]) 

title(’BANDSTRUCTURE OF GaAs/Al1GaAs (A1%25) SUPERLATTICE’ ) 
xlabel(’k (first BZ), 1/m'‘) 

ylabel(’Energy, eV’) 


hbhidbdoe upnomnpanoeendoonnaneqrennipeenene annie: 
LELEELIYLIRIGIISLIAIGSESLELEEESEELSEEEESESETEBEEVESESSEES 


EELELLIELELILLISERLESEEEEESSESELSS SESSBVYRBLIVEBLSSESSSSB 
ZESSEERESZ G. BASTARD’S DISPERSION RELATION S%%3%S%%S3%% 
ELTLLRETL ALL IFALSELEELELEESELELEELESEEBEEVEBSBELEEESEEEBSS 


rmmi=mi* (0.511e6) /9el6; 
mm2=(mi+0.0838*x)*(0.511e6) /9e16; 


wl=sqrt (2*mm1) /hbar; 
w2=saqrt (2*mm2) /hbar; 


dq=1000000; 
dE=0.000002; 


a=h+L; dz=d/600;D=h/2+L; 


it=ceil(s/dq); 
29-929 29S 99999 RSS R LEIS EEEE ESSE SEESESEESESSSESSEBESSESS 


30 








ESEESESSESES DISPERSION RELATION ESESSEESEBES 
a. SES ZLLRLERSESHS SEETEBEESESSEEESELESESESSESS 


oe 


n=1:1:200000; 
.05+n.*dE; 
“ae (cos ((sqrt(E)). *w1*L)) .* (cosh( (sqrt (Ec-E) ) .*w2*h)); 
p2=1/2* ((-sqrt((Ec-E) ./E) .*sqrt (m2/m1))+(1./ (sqrt ( (Ec- 
E)./E).*sqrt(m2/m1)))); 
p3=(sin((sqrt(E)).*w1l*L)) .* (sinh( (sqrt (Ec-E)) .*w2*h)); 


p=pl-p2.*p3; 


VILEEVLLEEVESEFSESVSSESSSSSSLSBESBEBLSVSSVBBVEBBSBSERB 
EEEZEEUSESESEEZSE FIRST BROULINE ZONE (BZ) %%83%3%%%3SE3%S 
ESSESLESESEESLESESSSSEESESSESEESTSSESES EEUSSSESE TEER CELE 


for v=1:it 
— gq=v*dq; 
g (v) =cos (q*d) ; 
[yy hh]=min(abs(g(v)-p)); 
Ee(v)=.05+hh*dk; 


aq(v)=q; 
end 
ff=(-1)*fliplr (qq); 
k=[ff£ qaqj; % FIRST BRUOILLOIN ZONE 
ee=fliplr (Ee); 
EE=[ee Ee]; % E(k) 


figure(2),plot(k,EE,’b.’) 
axis([-s s 0 0.45]) 
title([’BANDSTRUCTURE OF GaAs/AlGaAs (x=.25) SUPERLATTICE eee ]) 


xlabel(’k (first BZ), 1/m’) 
ylabel (‘Energy (eV)’) 
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APPENDIX B: COMPUTER CODE IN MATLAB PROGRAMMING 
LANGUAGE FOR NUMERICAL CALCULATION OF THE 
WANNIER FUNCTIONS | 


format long g 


EELELSLEVLSESEEEESEEEEEEEBSELEEEEEEEEESTEGESSEEEESEESEEEESESESEESEBEESS 
ESSESESSESSEEEESESSEEED INPUT PARAMETERS %&33 233333 24SS3SE3EES%SSS% 
ELELEELESLYELE SLE LSS EEASESEESEESEESEESEESEEEESEEESEEESLESESESEESESES ESSE 


EgG=input (‘Energy gap for GaAs (eV)’); 
hh=input(‘barier thickness (nm)'); 
LL=input (’well thickness (mm) ‘); 
x=input(’Al concentration (x)’); 


ELSELLELEILIE LISLE SSE EEELSESSSEEEESESSEBEESSESELEEGEESEEEESESSESEESESSES 
SELESEXLFSLESEALERSE DEFINING THE CONSTATNTS B8SSSESSSSSBVSSESVESESSF 
EEZLEESEERELEESELEEESEREELELELEELEELESEELETEEEKEGESESESTESSEEEESEEESS 


hbar=6 .5821e-16; 
ml=0.067; 
me=(0.511e6) /9el16; 
m2= (m1+0.0836*x) ; 


EgA=EgGt+1 .247*x; %AlGaAs bandgap for given (x) 
Ec=(0.65* (EgA-EgG) ) ; Conduction band offset 


h=hh/1000000000; 

L=LL/1000000000; 

a= (h) +L; 

D=(h/2)+L; DH=D+h; DHL=DH+L ;DHLH=DHL+(h); DHLHL=DHLH+L; 
s=3.14/d; 

ind=210; 

dz=a/70; 

ing=40; 

dgq=floor(s/ing) ; 

w=- (hbar%2) / (2*me* (dz%*2)); 


EESESEREERSESEEEEEEEEEEEEETEELEBSELEELESEETEEESTEEEELSEEETEESES ESSE SES SES 
EEESEELEESESESEELESSSS HAMILTONIAN {(H) EESESEESELSESESEEEGESS SES 
ESESLSESEESEESEESESEESESTEBTESESEEBESSESSESEESLSESESESSESTESESSESBESESSES 


for l=1:ind 
z(l)=d2z*1; 
L£ z(1)>(h/2) & z(1)<D 
V(1)=0; 
HH (1)=V(1)-(2/m1) *w; 
m(1)=mi; 
Ww (1)=w/ml1; 
elseif z(1)>D &z(1)<DH 
V(1)=Ec; 
HH (1) =V (1) -(2/m1) *w; 
Ww(l)=w/ml; . 
elseif z(1)>DH & z(1)<DHL 
V(1)=0; 
HH (1) =V(1)-(2/m1) *w; 
m(1)=m1; 
Ww (1)=w/ml1; 
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elseif z{(1)>DHL & z(1)<DHLH 
V(1)=Ec; 
HH (1) =V(1)-(2/m1) *w; 
Ww (1) =w/ml1; 
elseif z(1)>DHLH & 2(1)<DHLHL 
V(1)=0; 
HH (1) =V(1)-(2/m1) *w; 
m(1)=mi; 
Ww (1)=w/ml; 
else $z{(1})>DLHL %*& z{1})<DLHLD 
V (1) sEc; 
HH (1)=V(1)-(2/m1) *w; 
Ww (1) =w/ml1; 
. end 
end 


Ritw (floor (ind/18))=(2*w) / (ml+m2) ; 

Sww (floor (ind*7/18))}=(2*w) /(ml+m2) ; 
gww (floor (ind*11/18) )=(2*w) / (mi+m2) ; 
aww (floor (ind*13/18))={2*w) /(mi+m2) ; 
ww (floor (ind*i7/18))=(2*w) / (m1l+m2) ; 


HH1=diag (HH) ; 
W=Ww (1: (length (Ww)-1)); 
Wl=diag(W,1); 
W2=diag(W,-1); 


HH2= (HH1+W1+W2); 
nl=zeros (ind, 1) ;n2=zeros (ind,1);n3=zeros({ind,1); 


for k=1:2*ingq; 
if k>ing 
ql (k-inq) = (k-ing) *dq; 
Ql=(w/m1) * (cos (ql (k-ing) *d) -i*sin(ql (k-ingq) *d) ); 
Q2=(w/ml1) * (cos (ql (k-ing) *d) +i*sin (qi (k-ing) *d)); 
else 
ql (k) =-k*dq; 
Ql=(w/m1) * (cos (ql (k) *d) -i*sin(ql (k) *d)); 
Q2=(w/m1) * (cos (ql (k) *d) +i*sin(ql(k) *d)); 
end 
O11l=diag(Q1, (ind-1)) ;Q22=diag(Q2,-(ind-1)); 


H= (HH2+011+022) ; 

[jj kk])=eig(H); 
ee=sort (diag (kk) ); 

[£1 glj=find (ee (1) ==kk) ; 
[£2 g2]=find (ee (2) ==kk) ; 
[£3 g3]=find(ee(3)==kk) ; 


nl=n1+(jj(:,£1)); 
n2=n2+(jj(:,£2)); 
end 


EISESYYLLI LLL ILE LISLE ELLE LESLESEESEEEEESEESEBESSEESESESTESSESESSESESES 
TEXSEEAEYZELEEZ PLOTTING OF THE WANNIER FUNCTIONS %&3%S84SSSSEEESBSSS% 
EESEISEEEEESEEEBSBESEEEEVESYEEESEEETEEEESEESTEEESEEESEBTESEESESEESESSS 


g=1:ind; 

figure(1),plot((g.*dz) .*100000000,n1,’k’) 

title(’WANNIER FUNCTION OF GaAs/AlGaAs(A1%25) SUPERLATTICE’ ) 
xlabel(’‘z , nm ’) 
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ylabel(‘Amplitude (arbitrary) ‘) 
figure (2) ,plot((g.*dz) .*100000000,n2, ‘k’) 
title(’WANNIER FUNCTION OF GaAs/AlGaAs(A1%25) SUPERLATTICE ‘') 


xlabel(’z , nm ’) 
ylabel(’Amplitude (arbitrary)’)’ 
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